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Abstract 

Many structural transformations involve a group-nonsubgroup relationship between the initial 
and transformed phases, and hence are beyond the purview of conventional Landau theory. We 
utilize a systematic and robust methodology to describe such reconstructive martensitic transfor- 
mations by coupling group-theoretical arguments to first-principles calculations. In this context we 
(i) use a symmetry-based algorithm to enumerate transformation paths, (ii) evaluate the energy 
barriers along these transformation paths using all-electron first principles calculations, (iii) deduce 
the full set of primary and secondary order parameters for each path to establish the appropriate 
Ginzburg-Landau free-energy functionals, and (iv) for each path, identify special points of the pri- 
mary order parameter, as a function of local distortions, corresponding to the end product phase. 
We apply this method to the study of a pressure driven body-centered cubic (bcc) to hexagonal 
close-packed (hep) transformation in titanium. We find a generalization of the Burgers mechanism, 
and also find that there is no energy barrier to this transformation. In fact, surprisingly, we also 
find a region of volumes in which the intermediate path becomes more stable than either of the 
end-points (bcc or hep). We therefore predict a new orthorhombic phase for Ti between 51 and 62 
GPa. 

PACS numbers: 81.30.Kf, 64.70.Kb, 05.70.Fh 
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I. INTRODUCTION 



Polymorphism, i.e., the changing of a sohd from one crystal structure variant to another, 
is abundant in both nature and in engineering technologies. The body-centered-cubic (bcc) 
to hexagonal-close-packed (hep) structural transition discussed in this paper is a well-known 
example. For instance, iron undergoes this transformation at very high pressure (~ 13 GPa) 
in the core of the earth as well as under shock conditions in diamond anvil experiments^. 
Certain seismic activity has been attributed to strain generated during this phase change. 
Toughness and ductility in Ti and Zr, which undergo this transformation at high tempera- 
ture, are two important issues for the aerospace industry. In addition to these three elements, 
the bcc-hcp transformation is found in eighteen other elements, such as Ba, Mg, Hf, Sc, Tl, 
and Y, and their alloys. 

Polymorphic diffusionless transformations have been broadly classified as either displacive 
or reconstructive. In the displacive transformation there is no change in the first coordina- 
tion of the atoms in both the initial and the transformed phases. The energy differences 
between these two phases arise from changes in the secondary coordination. In terms of 
crystal symmetry, the transformed phase is a slightly altered derivative crystal structure, a 
subgroup, of the initial phase. In the reconstructive transformation, bonds are broken and 
then reformed to suit the transformed phase. Consequently, there is a change in the primary 
coordination. The transformed phase is no longer a slightly altered structure of the initial 
phase; it is not a subgroup structure of the parent phase (or, vice- versa). In diffusionless 
polymorphic transformations, individual atoms execute ordered and correlated motions, and 
significant transformational strains may also develop in this process. These are the so-called 
martensitic transformations (MT). 

The study of MT has a long and rich history. Early attempts by metallurgists to under- 
stand both displacive and reconstructive MT's ignored the need to maintain coherency be- 
tween the initial and transformed structures as the transformation progressed. That neglect 
hindered the development of an atomistic understanding of the transformation mechanism 
in terms of simple shears of the crystal cell coupled with local atom readjustments (shuffles). 
Thus, martensite theories often chose the phenomenological route of geometrically relating 
the initial and transformed phases, ignoring the atomistic pathways and microscopic mech- 
anisms driving the transformation. A well-known phenomenological theory of Wechsler, 
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Lieberman and Read successfully demonstrated that a MT can be understood at this level 
in terms of three basic deformationsiHS. Nevertheless, within the above mentioned phe- 
nomenological approach, microscopic mechanisms were identified for many reconstructive 
transformations, purely from the orientational relationships between the initial and trans- 
formed phases. In Table | we give examples of several common mechanisms. 

More recently, physicists have successfully used the phenomenological Landau theory to 
model displacive MT for a nearly continuous {i.e., weakly first order) transformation. Be- 
cause the displacive MT is nearly continuous, the invariant free-energy (invariant under the 
symmetry group of the initial phase) defines an order parameter (OP) onset. The Lan- 
dau philosophy has recently been extended to a phenomenological theory of reconstructive 
phase transitions by expressing the OP function as a density-wave function of the atomic 
displacement Note that an effective one- dimensional GmzhnTg-La.nda.u model for Ti and 
Zr§ has been combined previously with electronic structure calculation^. The bcc-hcp 
transition in titanium occurs around 1150 K. However, here we attempt to model this tran- 
sition at T = by varying the pressure: bcc is the high-pressure phase and the hep the 
low-pressure phase. The bcc-hcp transformation pathways and mechanisms elaborated be- 
low are generic to the bcc-hcp transition in a wide variety of elements and alloys. In effect, 
the temperature variation in titanium is mimicked by pressure variation at T = capturing 
the essential pathways of the transition, specifically the intermediate orthorhombic phase. 
As discussed below, we surprisingly find that one of the pathways becomes more stable than 
either the bcc or hep end-points in Ti for a range of volumes. We therefore predict a new 
orthorhombic phase transformation in Ti at low temperature. This paper is organized as 
follows. In the next section we introduce the concept of a common subgroup. In Sec. Ill we 
enumerate the various symmetry based pathways for the bcc-hcp phase transition. In Sec. 
IV we describe the first-principles calculations performed here. Section V gives an overview 
of the Landau free-energy fitting. Results and discussions are contained in Sec. VI. 



II. THE COMMON SUBGROUP PARADIGM AND TRANSFORMATION 
PATHS 

The bcc-hcp transformation in Ti considered here is a reconstructive martensitic trans- 
formation. There is no group-subgroup relationship between the space groups of the bcc 
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(Jm3m) and hep {PG^/mmc) phases partieipating in the transformation. BurgersQ origi- 
nally proposed the orientational relationship between the bcc and hep phases in Zr: the 
basal {OOOljhcp plane is parallel to a {llOjbcc plane, and a (1120) close-packed direction in 
the hep basal plane is parallel to a cubic (111) direction. A combination of three lattice 
distortions can lead to this transformationlilHEi: (i) a shuffle displacement of the {llOjbcc 
planes in the (110) directions. A displacement amplitude of ^ times the bcc lattice constant 
leads to the exact hep stacking sequence, (ii) A shear such as (112) [111] and (Il2)[lll] to 
squeeze the bcc octahedron into a regular hep one. The angle of the hep face in the basal 
plane changes from 70.53° to 60°. The Burgers mechanism was painstakingly determined 
by studying the orientational relationships between bcc and hep phases. It is understood in 
terms of the three parameters, shuffle, shear, and volume dilatation. 

In this paper we present a more general and robust approach for identifying reconstructive 
transition mechanisms based on crystallographic group theory. It is equally applicable to bcc- 
hcp MT's in other elements, or for other reconstructive MT's such as the transitions fcc-hcp, 
hcp-cu, etc. The distortion of a crystal can clearly destroy some of its symmetry elements and 
force it into one of its lower symmetry subgroups. As atoms move from one structure toward 
another, symmetry is lowered as the atoms leave their high symmetry locations. We assume 
initially that this lower symmetry is maintained until at some suitable, large displacement 
it is able to lock-in to the symmetry of the transformed phase. We thus hypothesize that 
the reconstructive transformations proceed via a subgroup common to both the initial and 
transformed phases. This requires the atoms to sensibly map from the initial structure 
onto the subgroup, and then onto the final structure. By imposing some restrictions on 
cell size changes the number of available transformational paths equals a finite number of 
maximal subgroups^, with compatible atomic mappings, common to both the initial and 
the transformed structures. We can then look for even lower symmetry structures along the 
path by testing stability within the maximal subgroup phase. The actual path chosen is 
ultimately decided by the nature of the interatomic interactions in the material. Coupling 
the symmetry with interatomic force calculations gives us a general, systematic, and robust 
recipe to explore the different pathways, and identify the minimum energy path(s) among 
them. This method is amenable to automation by computer and can quickly search for 
transformational paths between any two of the 230 crystallographic space-groups. 

We illustrate an application of our procedure by coupling the full potential linear 
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augmented-plane-wave (FLAPW) method to crystal symmetry information in order to enu- 
merate all possible paths for the bcc-hcp transformation in titanium. Based just on symme- 
try conditions, we initially obtain six possible mechanisms for the bcc-hcp transformation in 
titanium. Of these six, one is energetically favored and is similar to the Burgers mechanism. 
Our mechanism adds an effective interplanar variance to the three parameters used in the 
original Burgers mechanism. In this work, we extend the Landau philosophy to reconstruc- 
tive MT's by including local site symmetry information. Our procedure is distinct from that 
given in Toledano et al.0 in that we do not assume a transcendental stationary density-wave 
description but instead structure our approach on common subgroup symmetry aspects and 
a more conventional Landau position-dependent OP. However, we do use specific intracell 
correspondences of atoms in the two structures, which is not contained in the conventional 
Landau procedure, and look for the onset of new symmetry elements at the end product 
phase. These features are in common with the description given by Toledano et al. We 
fit a Landau free-energy functional over the multi-component OPs in this transition. The 
free-energy functional is fitted using the first principles energy values obtained for the most 
plausible transformation mechanism. 

The bcc-hcp transition in titanium occurs around 1150 K. Here we attempt to model this 
transition at T = by varying the pressure: bcc is the high-pressure phase and hep the 
low-pressure phase. In this sense our modeling cannot be considered as realistic in capturing 
the (temperature induced) transition in titanium. However, two new high-pressure phases of 
titanium have been observed recently, namely 7-Ti which is a distorted hep structure at 128 
and 5-Ti which is a distorted bcc structure at 140 GPa@. Interestingly, both of these 
structures are orthorhombic with the space group Cmcm and the lattice deformations can 
be interpreted within the framework of our generalized Burgers mechanism. We also predict 
a new orthorhombic phase transformation between 51 and 62 GPa at zero temperature. In 
this context our modeling is directly relevant to pressure induced transitions in titanium. 
In addition, the bcc-hcp transformation pathways and mechanisms elaborated below are 
representative of the bcc-hcp transition in a wide variety of elements (e.g. Zr, Hf) and 
alloys. In effect, the temperature variation in titanium is mimicked by pressure variation 
at T = capturing the essential pathways of the transition, specifically the intermediate 
orthorhombic phase. 
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III. CRYSTAL SYMMETRY BASED ALGORITHM TO DETERMINE THE 
MECHANISM FOR THE BCC-HCP TRANSFORMATION 



Methods based solely on orientational relationships between the initial and transformed 
phases force one to "guess" the transformational mechanism. Even symmetry based infor- 
mation coupled with atomic site correspondence does not unequivocally reveal the transfor- 
mational mechanism. Our approach here is to use "physics" insight (interatomic forces based 
on first-principles or some other suitable semi-empirical method) to interrogate the various 
common subgroup based symmetry pathways to determine plausible paths based on their 
relative energies. We expect that the minimum energy path, under the transformational 
conditions, will be favored, and we initially consider maximal symmetry pathways. Overall, 
such an approach reduces the systematic search of paths to a smaller number. In this pro- 
cess, we can also identify the primary and secondary order-parameters, the Ginzburg-Landau 
invariants, and declare the special OP lock-in values corresponding to the transformed phase. 

To describe a reconstructive phase transition, one must address two important questions: 
(1) How are the atoms mapped from one structure onto the other? (2) What distortional 
pathway do the atoms take between these two structures? The mapping question deals with 
how the two structures are related. The path question is more difficult. It deals with actual 
atomic displacements and strains that occur during the phase transition. In an earlier paper 
we presented a systematic procedure for obtaining possible mappings (paths) for a recon- 
structive phase transition in two dimensions (2D)0, namely the square lattice to triangular 
lattice transition, based essentially on symmetry considerations. In that treatment, it was 
assumed the pathway between the two structures proceeded by means of an intermediate 
common unstable structure (either an oblique lattice or a centered rectangular lattice) with 
definite space-group symmetry G. Note that the square— ^triangle transformation is a 2D 
analog of the bcc-hcp transformation. If we are considering a reconstructive phase transition 
from a structure with space-group symmetry Gi to a structure with space-group symme- 
try G2, then our description is actually a two-step process, Gi ^ G ^ G2, where G is a 
subgroup of Gi and also a subgroup of G2 so that each step, Gi ^ G and G ^ G2, is a 
transition with a group-subgroup relationship. 

The first question stated above is related to the mapping of atoms through a subgroup 
and was systematically addressed by a four step algorithmS We describe this algorithm 
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using the bcc-hcp transformation as an example. However, it is general and applies to other 
reconstructive transformations. The tables of Stokes and Hatch0 were used to obtain the 
needed information at each of the following four steps. 

Step-1: Find the subgroups common to both the initial (bcc) and transformed (hep) 
phases. To begin, it is convenient to look at common subgroups arising from distortions 
corresponding to k points of symmetry. We find seventeen subgroups which are subgroups 
of both the hep and bcc phase. 

Step-2: Subgroups arise from microscopic displacements and/or macroscopic strains. 
Strain cannot be the primary OP since the bcc-hcp transition involves primitive cell dou- 
bling. Shuffle is our primary OP and the induced representation formalism is used to deter- 
mine which irreducible-representations (IR's) allow microscopic atomic x, y, z displacements 
from each parent group into the subgroups. Only fifteen subgroups are consistent with this 
requirement. 

Step-3: Atoms are at the 2(a) and 2(c) Wyckoff positions in bcc and hep respectively. 
Subgroups that relate bcc and hep phases must have compatible Wyckoff positions. Each 
acceptable subgroup is connected to a specific IR and mode form (OP direction). Thus, there 
is a precise atom to atom identification throughout the transformational path. This ensures 
consistency as the Wyckoff position changes (or splits) as the transformation moves from bcc 
to hep. After applying this constraint, we are left with fourteen subgroups corresponding to 
fifty-eight possible mechanisms. 

Step-4: We determine a specific (correlated) displacement mechanism for each IR and 
subgroup. All possible mechanisms for the bcc-hcp transformation based on this algorithm 
are given in Table ||. Of these, we restrict our attention to the most simple mechanisms 
such as a maximum cell size change of two and a PI or P2 (one parameter) OP direction. 
As can be seen from Table II, these two restrictions yield three mechanisms; 63(P2,2), 
N2 63 (PI, 2), and 63(P1,2). This notation gives the IR, the identifying number of the 
subgroup, the order parameter direction, and the cell size change. The atomic displacements 
for these mechanisms are shown in Fig. [l|. Note that the space group symmetry of the 
common subgroup is the same and the atoms go into the Wyckoff c postions for each of the 
3 paths. 

The labeling of representations is that of Miller and Love@, the space group number is 
according to the International Tables^, and the OP direction is that of Stokes and HatchS. 
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Specifically, N and H refer to the high symmetry points in the bcc Brillouin zone, namely 
the face center and corner, respectively. Note that a phonon anomaly has been observed at 
the N-point in the phonon dispersion curves of Ti, Zr and rEMM. It is conceivable that 
the H-point mechanism may be related to the small anomaly seen at the H-point phonon 
dispersion of Mgil. The crux of the problem is to couple the pathway information to the 
relevant physics in order to identify which are the most energetically favorable paths. 

We generalized the above algorithmic procedure by allowing multiple primary order pa- 
rameter mechanisms and considering additional user input for some important parameters 
such as allowed strain, nearest-neighbor distances, and unit-cell size change. This computer 
program is called COMSUBS and filters possible subgroups through user defined conditions. 
The program was first used in the description of the rock-salt to CsCl structural change in 
sodium chlorideii. We applied the COMSUBS algorithm to the description of the transition 
of interest here, the hep to bcc transition in titanium. 

For the low-temperature hep structure, Gi = PQ^/mmc with lattice parameter ai = 
2.645 A, ci = 4.11 A. For the high-temperature bcc structure, G2 = ImSm with lattice 
parameter 02 = 2.91 A. We used the following criteria: 

(1) We considered only those subgroups where the length of the lattice generators is 
5.03 A or less. This condition effectively limits the size of the unit cell in the common 
subgroup, which in turn limits the size of the allowed unit cell parameters. 

(2) We considered only those subgroups where the principal elements of the strain tensor 
are less than 1 + e and greater than (1 + e)^^, where e = 0.5. This condition limits the strain 
(lattice vector lengths) to be within this allowance and also the angle relationships between 
the group and its subgroup. 

(3) The nearest-neighbor distance is 2.56 A in Gi and 2.52 A in G2- We considered only 
subgroups where the nearest-neighbor distance in the structure halfway between Gi and G2 
is greater than 2.03 A (80% of the average of 2.56 A and 2.52 A). If the nearest neighbor 
distance is less than our chosen distance of 2.03 A we reject the path. Along a rejected path 
the lattice would need to expand to allow atoms to pass one another and we assume such a 
strain implies an energy barrier unfavorable to this path. 

(4) We considered only maximal subgroups. These define the possible mappings of atoms 
in Gi onto atoms in 6*2. Subgroups of these maximal subgroups do not introduce new 
mappings. They only alter the path by allowing additional distortions in G along the path. 
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We check for these additonal distortions when a maximal subgroup is selected. 

The criteria given above were chosen somewhat subjectively. Each of the criteria could 
be relaxed and the list would thus be extended. We will see below that the values we chose 
resulted a list which contained a subgroup that upon initial screening and additional checks 
for stability is the likely candidate for the pathway mechanism. 

Using these criteria, we obtained six subgroups from COMSUBS, which we list in Table 



Il|. For the first two entries in the table there is no change in the size of the primitive 



unit cell relative to the hep cell. The next two entries are subgroups where the size of the 
primitive unit cell is doubled. The last two entries are subgroups where the size of the 
primitive unit cell is three times larger than that of the hep cell. Subgroups with larger 
primitive unit cells were not found since we limited the length of the generators according 
to criterion (1) above. We label the six potential mechanisms we obtained as, (a) ort63A, 
(b)ort63B, two pathways via orthorhombic space-group 63, (c) monl4, a pathway via the 
monoclinic space-group 14, (d) monl5, via the monoclinic space-group 15, (e) mon9, via the 
monoclinic space-group 9, and (f) tri2, via the triclinic space-group 2. 

The list obtained by COMSUBS does contain groups in common with those obtained 
from the four-step algorithm described earlier (for example ort63B), some are excluded by 
COMSUBS since they do not meet the criteria listed above, and new pathways are obtained 
due to the allowance of coupled parameters (for example monl4 couples and A'"2"). 

Prescreening calculations showed (see Table [TTTI ) that the second pathway based on the 
orthorhombic space-group Cmcm(63) had the lowest energy barrier among all these paths, 
and it was chosen for detailed study. The calculations were done at midpoint structures as 
determined by the bcc-hcp endpoint structures. This distortional mechanism is described in 
terms of an appropriate high-symmetry point (N) of the bcc Brillouin zone and irreducible 
representation Nj, with alternating shuffles in the ±[011] directions in the (Oil) planes. 
This mechanism is closely related to the original Burgers mechanism. 



IV. MECHANISMS FOR THE BCC-HCP TRANSFORMATION IN TITANIUM 

The most favored transformation path should have the lowest energy among all the 
paths considered here. Energy of the different paths were computed using an all-electron, 
full potential linear augmented-plane-wave (FLAPW) method^. Calculations were per- 
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formed scalar-relativistically, neglecting spin-orbit coupling for the valence electrons. Local 
orbitalS were added to enhance the variational freedom and allow the semi-core 3s and 3p 
orbitals to be treated along with the valence electrons. The added energy parameter was 
used to simultaneously treat the residual s and p character of the valence electrons. Our 
results are insensitive to small changes in all energy parameters. 

Our calculations used the generalized gradient approximation (GGA) of Perdew, Burke, 
and Ernzerholii for the exchange-correlation functional. A muffin-tin sphere with a radius 
Rmt was used to define the augmented plane-wave basis functions. No restriction is made as 
to the shape of the potential or charge density. The size of the FLAPW basis was determined 
by a plane-wave cutoff, K^ax- The K^ax value is given by the relation RuxKmax = 9.0, and 
was found to be adequate in determining the system energy. 405 irreducible points in the 
first Brillouin zone of orthorhombic Ti were sampled according to the improved tetrahedron 
methodil for Brillouin zone (BZ) integrations. Total energies changed by less than 0.014 
eV/atom on increasing the number of BZ points sampled. All self-consistent calculations 
were iterated until the total energy changed by less than 0.1 meV/atom. 

Based on the FLAPW calculations, we found that an intermediate orthorhombic pathway 
becomes more stable than either of the bcc or hep end-points for volumes between 49.197 
and 50.235 A . The intermediate orthorhombic phase gives rise to orthorhombic-bcc and 
hcp-orthorhombic transitions at 61.8 and 51.3 GPa, respectively (See Fig. 0). Because the 
orthorhombic unit cell and the atoms within it were relaxed to find the lowest energy, these 



parameters change slightly in this range. In Table |^ we list the relevant atomic and unit 
cell coordinates near the middle of the stable orthorhombic region. The y coordinate refers 
to the value of y in the specification of the Wyckoff c position of the form 0,?/, 1/4. If 
we ignore this orthorhombic intermediate phase, we find that the bcc-hcp transition occurs 
at a constant pressure of ~ 57 GPa. The transition pressures were determined by the 
slope of the common tangent to the energy versus volume curves of the two appropriate 
phases participating in the transformation. Since the curves are very close to each other 
in energy, we are unable to determine from the common tangent construction whether the 
orthorhombic phase is more stable than a mixture of bcc and hep crystals. We do see that 
there are regions wherein the orthorhombic system has lower energy than either the hep or 
bcc phases. 

Since our preliminary calculations pointed to a mechanism operating via a common or- 
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thorhombic subgroup, we chose to investigate it in more detail. The orthorhombic cells 
were fully relaxed with respect to its lattice parameters and internal coordinates, a total of 
four variational parameters, to determine minimum energy configurations. Such minimum 
energy surfaces for the common orthorhombic subgroup were determined at five different 
volumes, 49.316, 49.434, 49.716, 50.054, and 50.235 A , in the transformation region. We 
find that the energy of the orthorhombic phase is also sensitive to its lattice constant along 
the c-axis. All of the five volumes possessed minimum energy when the orthorhombic cell 
had a length between 4.04 and 4.07 A. This c-axis value is significantly different from the 
value of 4.11 A corresponding to the clamped orthorhombic structure of the hep phase. 
Such a variation in the c-axis lattice constant was not considered in the original Burgers 
mechanism. The energy versus volume plot for the bcc, hep and orthorhombic phases in 
the vicinity of the transformation region is given in Fig. ^. We can see that there is no 
effective energy barrier to the bcc-hcp transformation when the transformation proceeds via 
the orthorhombic phase. 



V. DETERMINING FREE-ENERGY FUNCTIONAL FROM FIRST- 
PRINCIPLES DATA 

Here we demonstrate that it is possible to rigorously fit the Ginzburg-Landau free-energy 
functional to first-principles energy data. These functionals will be used in our future work. 
One primary and three secondary order parameters (OP's) completely describe the bcc-hcp 
transition through the orthorhombic phase. They are the intracell atom "shuffle" (?7i), the 
shear {i]2), deviatoric (173), and volumetric (774) strains. Shuffle is the primary OP and is in 
units of Bohr. The orthorhombic c-axis parameter gives us a measure of ri2, and is also in 
the units of Bohr. The deviatoric OP rj^, defined as the ratio of a- to 6- axis parameters of 
the orthorhombic cell, is dimensionless. The volumetric OP (r/4) has the units of Bohr'^. The 
full form of the Ginzburg-Landau free-energy functional, involving all four OP's, is given by 
the equations: 

FLivi, V2, Vs, V4) = Mvi' + Mvi^ + -43771*^ 
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+^137/4^ + Aiiri2^ri3 + Air^ris^ 
+AiQrii^rii + Ai'jri2^rii + An^rj^^rji, 



(1) 



Fcirii) = giVi/ + 92{Vi/ + Vi/) + 93Vi,yVi,z- (2) 

Equation (P is the Landau portion of the free energy and gives the energy for a homoge- 
neous phase arising from coupled primary and secondary OP contributions. Equation 
contains the (nonlocal gradient) Ginzburg portion of the free energy and gives the energy of 
heterogeneous solutions, e.g., energies of domain walls configurations, etc. The notation in 
Equation (0) indicates derivatives of the OP components, e.g., rji^^ is the x derivative of rji. 

It is important to note that this free-energy is at T = K. The coefficients of Fi, the 
Landau portion of the free energy given above, were determined by a Powell non-linear least 
square fittingil of the first-principles energy data, along the bcc-hcp transition path. Since 
the purpose of this work is to demonstrate the feasibility of fitting GL free-energy functional 
to first-principles energy data, we carried out the fit at one fixed volume (50.054 A^) falling 
within the transformational volume region and using 58 energy values corresponding to this 
volume obtained from first-principles calculation. The coefficients are given in Table |V]. The 
fitted energy values had an error of less than 0.5 %. Here, the summation index i is over all 
the 58 energy values. 



VI. RESULTS, DISCUSSION AND CONCLUSION 

As introduced earlier, Burgersi originally derived the orientational relationship between 
the bcc and hep phases for the bcc-hcp transformation. This orientational relationship was 
translated into a mechanisnJiiBlii with three parameters: (a) shuffle of the atoms in the 
{llOjbcc planes in the (110) directions; a displacement amplitude of ^ times the bcc lat- 
tice constant leads to the exact hep stacking sequence, (b) a shear such as (112) [111] and 
(ll2)[lll] to squeeze the bcc octahedron to transform it into a regular hep one; this changes 
the angle of the hep face in the basal plane from 70.53° to 60°, and (c) an appropriate volume 
dilatation to reach the correct end volume. Such a mechanism is similar to our pathway via 
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the intermediate orthorhombic subgroup in the above three respects. However, our pathway 
has an additional characteristic in that it also allows the orthorhombic lattice constant along 
the c-axis to vary. In retrospect, although such a variation is intuitive, it is surprising that 
Burgers apparently did not factor such a possibility into his mechanism. We suggest that the 
material should follow our more general four parameter pathway during the transformation. 
In essence, the utility of our coupled symmetry and first-principles calculations is demon- 
strated in that we find a Burgers like mechanism, albeit a more generalized one. We used 
a general and robust algorithm to ferret out the correct mechanism, one that is consistent 
with the symmetry and physics of the transformation. There is no energy barrier for this 
new generalized Burgers mechanism. 

We have carried out our calculations for T = which may be more accurate for pres- 
sure induced transitions such as hep —^ooin titanium. However, bcc — > hep is a thermal 
transformation which takes place around 1150 K. Therefore, realistically finite temperature 
phonon effectsU cannot be overlooked. We, however, calculated the energies of various 
frozen phonon structures (using the computer program FROZSL-INIT by L.L. Boyer and 
H. T. Stokes) that deviate from our calculated orthorhombic pathway and list the relative 
energies of these distorted structures with respect to the reference orthorhombic energy in 
Table |VI[ It is clear that the energy of these phonon distorted structures are higher than 
that of our orthorhombic pathway. Thus, we expect the overall transformation path to 
remain Burgers-like at finite temperature except that energy barriers for competing paths 
may change (~100 meV for this temperature). Work is underway to add the effects of 
entropy in the free-energy. In addition, defects are invariably present during many phase 
transformations. Depending on the defect energetics and concentration they may alter the 
transformation path or in some cases arrest it altogether (e.g. a few ppm oxygen arrests the 
shock-induced hep uo transformation in titanium). For low energy and low concentration 
of defects we do not expect our results to change significantly. However, there is a need for 
systematic evaluation of how defects affect transformation paths and barriers. 

To our knowledge this is the first systematic attempt to enumerate transformation paths 
based on symmetry and then to integrate this information with first principles calculations 
for the bcc to hep transformation. Our approach is general and mathematically rigorous 
with no ad hoc or empirical elements involved. We believe that the unique combination 
of symmetry and electronic structure is an equally powerful tool to identify transforma- 
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tion mechanisms in many other reconstructive transformations in nature, e.g., the Bain 
mechanism for fee to bee transformation!, the Wentzeoviteh-Lam mechanism for fee to hep 
transformation!, and the diamond to NaCl structure. As we have seen for Ti, this ap- 
proach also has the possiblility of predicting new, unforseen phases that would be difficult 
to find by ad-hoc guesses. Finally, in addition to a Burgers-like mechanism we predicted 
new mechanisms, N2 and H4, which may well exist in other elements and alloys including 
Mg. 
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FIG. 1: Three simplest mechanisms that take a bcc phase into a hep structure, and based on a 
common orthorhombic subgroup (Space-Group 63, Cmcm): (a) NJ, (b) HJ, and (c) • 
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FIG. 2: Energy versus cell volume for the bcc, hep and orthorhombic phases of titanium. The 
orthorhombic phase has the lowest energy within the transformation volume region. 



TABLE I: Reconstructive transformation mechanisms from the literature. 



Transformation 


Mechanism 


bcc-fcc 


Baini 


bcc-hcp 


BurgerJi 


bcc-w 


Cook§ 


fcc-hcp 


Wentzcovitch-Lamll 


hcp-w 


SilcocJ 
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TABLE II: List of all 58 possible mechanisms for the bcc-hcp transformation allowed by symmetry. 
For each IR, we specify the common subgroup, OP direction, and cell size change. The notation 
for order parameter direction is that of Stokes and Hatchlli. 



IR 


Subgroup (OP direction, Cell Size Change) 




11 (Cl,2), 63 (P2,2) 




i (o4,4j, i ^DJJi,oj, Z [(^Z,^), Z [DO,o), z^4JJi, oj 

5 (S12,8), 5 (4D3,8), 8 (4D6,8), 9 (S13,8), 9(4D4, 8) 
12 (C8,8), 15 (S6,8), 42 (S14,8), 43 (C19,8), 
59 (Cl,4), 63 (Pl,2) 


AT~ 


1 tQ'i A\ 1 (RT^^ fi\ 9 (no a\ 9 /Q/i 9 /'/im 

i (^00,4j, i ^DUi,oj, / (^V^Z,4J, Z (^o4,oj, Z (^4JJi,oJ 

5 (C8,4), 5 (S12,8), 5 (4D3,8), 5 (4D6,8), 8 (S13,8) 
8 (4D4,8), 12 (C9,8), 12 (S6,8), 42 (C19,8) 




1 (S4,4), 1 (6D1,8), 2 (C2,4), 2 (S3,8), 2 (4D1,8) 
5 (S12,8), 5 (4D3,8), 8 (S13,8), 8 (4D4,8) 9 (4D6,8) 
12 (S6,8), 15 (C8,8), 42 (C19,8), 43 (S14,8), 
58 (Cl,4), 63 (Pl,2) 


P4 


1 (6D,4), 2 (S3,4), 5 (S12,4), 5 (S18,4), 5 (4D1,4) 
8 (4D6,4), 9 (S5,4), 12 (C8,4), 15 (C2,4), 44 (C22,4) 



TABLE III: Energies from initial pre-screening calculations for six mechanisms given by the COM- 
SUBS algorithm. Calculations were done at midpoint structures as determined by the bcc-hcp 
endpoint structures. Energy is in eV/atom, and the lattice-parameters are in A. ^ denotes an 
unrelaxed energy value. All of the rest are relaxed energy values. 



Mehanism 


Lattice-Constants 


Cell- Angles 


Energy 


ort63A Cmcm 


(2.78, 4.35, 4.11) 


(90, 90, 90) 


-4.51 


ort63B Cmcm 


(2.78, 4.35, 4.11) 


(90, 90, 90) 


-5.05 


monl4 P21/c 


(2.57, 4.33, 4.95) 


(90, 116.1, 90) 


-4.84 


tri2 PI 


(4.77, 4.63, 2.54) 


(85, 87.2, 62.5) 


-4.43 


mon9 Cc 


(10.25, 4.32, 4.47) 


(90, 48.7, 90) 


-4.70t 


monl5 C2/c 


(8.18, 4.27, 5.90) 


(90, 46.2, 90) 


-4.24 
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TABLE IV: Sample atomic and unit cell coordinates for the intermediate orthorhombic phase. The 
unit cell parameters are in A, the cell angles are in degrees, and the y coordinate is in fractional 
atomic units. The y coordinate refers to the value of y in the specification of the Wyckoff c position 
of the form 0, y, 1/ 4. 



Space-Group 


Lattice-Constants 


Cell- Angles 


y coordinate 


63 Cmcm 


(2.92, 4.29, 4.02) 


(90, 90, 90) 


0.295 



TABLE V: Ginzburg-Landau free-enegy coefficients in Eq. 1 for the bcc-hcp transformation in Ti 
via the N^" mechanism. Fl is in eV, rji and 772 are in A, 773 is dimensionless, and 7/4 is in A^. 





= -48106.70 


A2 


= +319.26 


^3 


= +1412.71 


Aa 


= -312.31 


As 


= -848.33 


Ae 


= +6.94 


A7 


= +26707.92 


As 


= +7434.47 


A9 


= -44094.63 


Aw 


= -19221.73 


An 


= -2145.13 


Au 


= +0.00 


Ai3 


= +0.00 


Au 


= -565.07 


Ai5 


= +27635.83 


Aw 


= +0.00 


An 


= +0.00 


Ais 


= +0.00 
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TABLE VI: Relative energy (AE), in units of meV/atom, of various frozen-phonon structures that 
deviate from the transformation pathway via the common orthorhombic subgroup. The energy of 
the orthorhombic subgroup is used as reference. The cell volume is in A^/atom. Some IR's have 
multiple frozen phonon configurations. Energies of IR's with multiple frozen phonon configurations 
are listed in the same line. 



IR 


Snarp-PTonD 


AE 


r+ 
^ 1 


R '\ (! ID ^TTI 
\JfJ 1 1 Vl^l I 


n 




11 P121/m 


170 


3 


12 C2/m 


1775 


Y+ 
^ 1 


..J _L -X. r 1 V 1 t VKv 


792 




± loll vUi 


9219 


-^3 


IlS Pt) J) TY) 
(JU J. 1 LI vl 1 1/ 


71 


z 


62 Pnma 


639 




57 Pmmm 


391 


Y4 


59 Pmmn 


31 


Zi 


12 C2/m 


41, 449, 255 


Z2 


15 C2/c 


54 


St 


11 P121/m 


187, 551, 650, 305 


ST 


14 P21/c 


332 


S2 


11 P121/m 


347, 483, 612 


Ti 


12 C2/m 


738, 622, 364, 238 
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